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Новые компьютерно-ориентированные формы 
математических моделей авиационных 
газотурбинных двигателей 


Для круга задач оптимального управления, оценки технического состояния и диагностики 
авиационных газотурбинных двигателей рассмотрены вопросы усовершенствования известных 
кусочно-линейных дифференциальных и построения новых интегральных математических моделей, 
адекватно отражающих статику и динамику объектов управления и диагностики. Предложены новые 
математические модели в виде интегральных уравнений Вольтерра и методы их компьютерной 
реализации. 


Введение 


Общая постановка проблемы построения адекватных имитационных и диагности- 
ческих моделей статики и динамики авиационных газотурбинных двигателей (ГТД) 
обусловлена как важностью теоретического исследования совокупности термогазоди- 
намических процессов в ГТД, так и практической значимостью ожидаемых результатов 
для решения задач оптимального и адаптивного управления, диагностики и оценки 
технического состояния в течение жизненного цикла объекта. Газотурбинный дви- 
гатель является одним из наиболее совершенных и наиболее сложных объектов 
энергетики, в котором рационально объединено взаимодействие термогазодинамических 
процессов для достижения высоких технико-экономических и экологических пока- 
зателей. Жизненный цикл ГТД включает этапы проектирования, создания опытных 
образцов, испытаний и эксплуатации. Информационным обеспечением всех этапов 
жизненного цикла ГТД является его математическое моделирование. Компьютерно 
реализованные математические модели (ММ) высокого уровня сопровождают 
двигатель в течение его жизненного цикла, представляя собой его информационный 
портрет. Поэтому проблема усовершенствования математических моделей в направ- 
лении повышения уровня адекватности с допустимым уровнем сложности имеет 
важное научное и прикладное значение. 

Следует отметить, что применение математических моделей на этапе разработ- 
ки в рассматриваемой отрасли передовых технологий вошло в практику отработки 
новых изделий, позволяет существенно сократить объемы стендовых испытаний, 
заменяя ряд их дорогостоящих этапов компьютерным моделированием [1]. Соответст- 
венно сокращается продолжительность этапа испытаний и число опытных экземпляров 
изделий, что обеспечивает значительный технико-экономический эффект от приме- 
нения ММ. С другой стороны, в процессе эксплуатации ГТД необходимость применения 
ММ для замещения сигналов измерительных каналов при их отказах уже введена в 
нормативные документы [2]. В работах авторов [3], [4] и других работах, обосно- 
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вано предложение создания электронного паспорта ГТД в виде его индивидуальной 
ММ, которая сопровождает двигатель в течение его жизненного цикла, обеспечивая 
возможность эксплуатации по техническому состоянию. Поэтому проблема усовер- 
шенствования и построения новых ММ для авиационных ГТД имеет теоретическую 
значимость и прикладную востребованность. 


1. Постановка проблемы и цель исследования 


Применяемые ММ авиационных ГТД можно разделить на термогазодинами- 
ческие, поузловые, поэлементные и кусочно-линейные динамические [5]. Детальная 
классификация и анализ применяемых ММ ГТД приведены в [6]. 

Среди моделей авиационных ГТД самостоятельное значение и существенную 
важность имеют имитационные ММ статики и динамики, для которых ГТД пред- 
ставляется единым объектом с совокупностью внешних воздействий и выходных 
координат. Имитационные модели ГТД в настоящее время являются необходимой 
составной частью ТЗ на ЭСУ профиля ЕАПЕС и используются в стендах-имитаторах, 
позволяя на этапе проектирования совершенствовать алгоритмы контроля состояния, 
управления режимами и диагностики [3], [4]. Как правило, имитационные ММ представ- 
ляют собой совокупность нелинейной статической и кусочно-линейной динамической 
моделей (КЛДМ) [7], [8], таким образом, в каждой точке рабочего режима ММ является 
линейной. Успешность применения КЛДМ обусловила разработку ряда их модификаций, 
среди которых наиболее важными представляются стохастические Марковские ММ [7] 
и ММ в виде векторно-матричных дифференциальных уравнений с неопределен- 
ными собственными значениями [9]. 

Рассматривая круг нерешенных задач, препятствующих широкому внедрению 
ММ на всех этапах жизненного цикла ГТД, следует отметить важнейшие из них, а 
именно: 

— имеющиеся совершенные термогазодинамические ММ непригодны для решения 
задач адаптивного и оптимального управления режимами ГТД и пока не могут быть 
реализованы в режиме реального времени в электронных системах управления; 

— применяемые дифференциальные ММ ГТД имеют феноменологический характер, 
не отражают физику процессов в объекте и плохо приспособлены к компьютерной 
реализации в бортовой аппаратуре. 

Таким образом, известные модификации ММ имеют недопустимо высокий 
уровень сложности при компьютерной реализации, что затрудняет их прикладное 
применение. Практика применения ММ ГТД требует отыскания таких их эквива- 
лентных форм, которые бы без снижения уровня адекватности имели допустимую 
компьютерную сложность, имея в виду в перспективе реализацию в бортовой ЭСУ. 
Разработка новых ММ ГТД, имеющих более общий характер и известные [10] пре- 
имущества при численной реализации, позволяет повысить уровень адекватности 
применяемых ММ и приблизить уровень их алгоритмической сложности к возмож- 
ностям бортовых средств. 

Целью настоящей работы является усовершенствование известных кусочно- 
линейных динамических дифференциальных моделей на основе их модальной фор- 
мы и обоснование новых ММ установившихся и переходных режимов ГТД в виде 
интегральных уравнений Вольтерра. 
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2. Основные результаты исследования 


2.1 Модальная форма математических моделей ГТД 


Переход к модальной форме ММ ГТД путем диагонализации векторно-мат- 
ричных уравнений динамики предложен Р.Л. Лейбовым [9] для решения задачи на 
неопределенные собственные значения сопровождающей матрицы характеристи- 
ческого полинома, которая сводится к совокупности задач линейного и нелинейного 
программирования. Размерность последней задачи столь велика, что теряет содержа- 
тельное значение для практики. Тем не менее, предложение перехода к модальной 
форме имеет весьма важное значение ввиду следующих ее преимуществ: 

— упрощения численной реализации и необходимых вычислительных ресурсов в п раз, 
где п — порядок системы, 

— достижения возможности качественного исследования ММ, ее фундаментальных 
свойств, а именно устойчивости, управляемости и наблюдаемости. 

Такой переход, в отличие от [9], предлагается осуществить таким образом, 
чтобы исключить весьма существенный недостаток КЛДМ - разрывность ускорений 
в точках сопряжения диапазонов, являющийся неблагоприятным фактором, особен- 
но при моделировании турбовальных ГТД. 

Для классической формы уравнений состояния, которыми описывается дина- 
мика ГТД в каждой точке рабочего режима, 


Е о а 
—х= Ах+ Ви, 
[и 
т и (1) 
у=Сх + ри, 
модальной является следующая их форма 
НЕ ы: 
—У=5у + Вии, 
[1 
(2) 


у=Су+ ри, 
где 5 =ГАГ`! — диагональная матрица из собственных значений 5, матрицы А, 
Г — матрица диагонализации Г = [Я т | ‚ являющаяся решением задачи на 
собственные значения 
А =. (3) 
Собственные значения 5, удовлетворяют уравнению 
4е((5Е — А)=0. (4) 
Поскольку обычно матрицы в (1) заданы некоторой совокупностью в зависимости от 


значений режимного параметра 7’, где ] - номер участка дискретизации, возникает 


естественное предложение их непрерывной аппроксимации. Теоретические основы 
такого представления известны в матричной алгебре в виде теории А -матриц [11], то 
есть матриц, элементы которых являются многочленами от параметра ХА. 

Можно видеть, что для предлагаемой аппроксимации имеются глубокие основа- 
ния. Действительно, трудно полагать, чтобы в реальной динамической системе при 
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изменении режимного параметра собственные значения перемещались по разрывным 
траекториям (случай переключения КПВ рассматривается отдельно). Такие траектории 
должны быть непрерывными и иметь точки бифуркации, в которых происходит измене- 
ние кратности и разделение на комплексно-сопряженные пары. В КЛДМ (1) координаты 
состояния представляют собой отклонения оборотов турбин от их статических 
значений в каждой точке рабочего режима, поэтому переход между соседними участ- 
ками аппроксимации сопровождается возникновением разрыва первого рода для 
производной оборотов, то есть ускорения, что противоречит физической сущности 
динамики ГТД. Поэтому форма (2), представленная через А -матрицы, является более 
совершенной аппроксимацией реальной динамики нелинейной системы и ГТД, в 
частности. 


Следуя теории А -матриц [4], если А(А) является регулярной матрицей 
(4её(4(1)) = 0), то она эквивалентна диагональной матрице с элементами, кото- 


рые являются инвариантными многочленами А(Л) . Эта диагональная д -матрица 


называется канонической формой Смита. Известный путь приведения А-матрицы к 
канонической форме Смита состоит в построении унимодальных матриц преобра- 
зования, которые, в свою очередь, представляются в виде произведения конечного 
числа элементарных Д -— матриц [9]. Громоздкость и неудобство для компьютерной 
реализации известного алгоритма затрудняют его практическое применение. Поэтому 
предлагается осуществить переход к модальной форме уравнений динамики ГТД на 
основе решения траекторной задачи на собственные значения в виде 


4(1)%, (1) =5% (1), (4), (5) 
которая эквивалентна следующей задаче 
А(Л)У (4) =7(1)5 (4), (6) 


где А(Л) является известной, 5 (1) определяется (4), а неизвестными являются траек- 


тории элементов собственных векторов. 
Если найдется решение траекторной задачи на собственные значения (6), то 
модальная форма уравнений динамики ГТД приобретает следующий вид 


а ыы =: 
2. = 5(1) + Ва (ли, 
ГА 
У=С,(4)7+Б(%)я, (1) 
а изменения координат СОСТОЯНИЯ определяются соотношением 


= ' (4). (8) 


Преимущества формы (7) достаточно очевидны, поскольку достаточно прозрач- 
ными становятся запасы устойчивости, в численном виде могут быть определены 
показатели управляемости и наблюдаемости в виде параметров обусловленности 


матриц Ву (4) и Су (1) на режимах, имеет место непрерывность ускорений. 


Необходимо отметить, что исследование траекторий собственных значений имеет 
самостоятельное значение, поскольку точки бифуркации таких траекторий, соот- 
ветствующие возникновению пары комплексно-сопряженных собственных значений, 
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позволяют выполнить качественный анализ ММ. Действительно, такая бифуркация 
означает потенциальную возможность возникновения колебательных режимов в газо- 
воздушном тракте и позволяет локализовать подсистему таких режимов, например, в 
виде: турбина 1 — турбина 2, или турбина 2 — турбина 3, что является важным для 
оценки запасов ГДУ и причин их изменения. 

Траекторная задача на собственные значения (6) имеет нетривиальный харак- 
тер и недостаточно полно рассмотрена в теории А -матриц. Первая проблема в виде 
комплексного характера матрицы диагонализации решается использованием блочно- 
диагональной формы Шура. 

Вторая проблема состоит в определении такой совокупности собственных векто- 
ров, которые бы образовывали непрерывные траектории. Стандартное программное 
обеспечение универсальных программных пакетов такой возможности не предо- 
ставляет. Однако такие проблемы носят непринципиальный технический характер и 
вполне преодолимы. Тестирование предлагаемой ММ выполнено на данных полетной 
регистрации трехвального ТРДД в режимах «приемистость — встречная приемистость» 
и его результаты иллюстрирует рис. 1. На рис. 1 представлены изменения оборотов 
турбин и их модельные значения по форме (1) С.В. Епифанова и определенные по 
модальной форме. Как это следует из результатов компьютерного моделирования в 
сопоставлении с реальными данными, различия не являются существенными. 
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Рисунок 1 — Результаты компьютерной реализации модальной формы в 
сопоставлении с реальными данными 
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2.2 Интегральная форма математических моделей ГТД 


Рассматривая решение (1) 
1 
4 (= [Ф(1,5)8(5) 48 (5) 45. (9) 
1, 


где Ф -— переходная матрица (1), можно заметить, что решение задачи на собст- 
венные значения и построение их траекторий в зависимости от режимного параметра 
по сути является аппроксимацией ядра интегрального преобразования (9) разност- 
ным экспонентным ядром, так как для каждого из участков аппроксимации матрицы 
А, В, С, О являются фиксированными 

Отсюда следует предположение, что основные термогазодинамические пара- 
метры ГТД могут быть связаны более общей, в сравнении с (1), (7) математической 
моделью, а именно ММ в виде интегральных уравнений Вольтерра. Такая ММ 
первоначально построена для установившихся режимов ГТД, а в последующем и для 
динамических режимов. 

Исходной гипотезой для построения ММ для установившихся режимов ГТД 
является предположение, согласно которому энергетический параметр двигателя в 
виде степени повышения давления, определяющий тяговые характеристики ГТД, 
аддитивно связан с оборотами его турбин 


М№ 

1=1 
а градиенты оборотов турбин по расходу топлива, в свою очередь, определяются этим 
энергетическим параметром 


4п; (Ст) _ 
АС’ -В.(бг я, (Ст), (11) 
1=1,М, 


где {а ь (х )} : } р, | (х)} — системы некоторых линейно независимых функций. 


Из (10), (11) непосредственно вытекает эквивалентная ММ в виде интегрального 
уравнения Вольтерра П-го рода с сепарабельным ядром 


=.(6,).а (бу) (,)* 


{ (12) 
У) (6) [В (5)л, ($)4 


н 

Идентификация ММ (12) выполнена по статическим характеристикам и базам 
данных испытаний того же двигателя и дала приемлемые результаты [5]. На рис. 2 
представлены нормированные значения изменения параметров двигателя и их модель- 
ные значения, определенные по интегральной форме. 

Среднеквадратическое отклонение по основным параметрам ГТД и, лк не 
превышает 1 %. Преимущества ММ (12) представляются достаточно значительными: 
основные характеристики ГТД определяются небольшой совокупностью линейно 
независимых функций, которые имеют ясное физическое содержание и довольно 
несложно определяемые в численном виде по экспериментальным данным. 
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Рисунок 2 — Результаты компьютерной реализации интегральной формы 
в сопоставлении с реальными данными 


Заключение 


Таким образом, практика проектирования современных электронных систем изме- 
рения параметров, управление и диагностики требует создания адекватных компьютерно- 
ориентированных математических моделей авиационных двигателей. Адекватность 
модели ГТД должна быть обеспечена ее соответствием реальным данным испытаний, 
а возможности компьютерной реализации — удовлетворять требованиям к бортовой и 
наземной диагностирующей аппаратуре. Предлагаемые новые ММ авиационных 
ГТД имеют ряд важных преимуществ в сравнении с известными моделями благодаря 
возможностям их численной реализации. Обучение разработанных ММ выполнено с 
применением современных алгоритмических и программных средств и позволяет 
индивидуализировать их параметры относительно конкретного экземпляра двигателя. 

Результаты исследований по созданию и усовершенствованию ММ ГТД, методов 
их компьютерной реализации позволили решить ряд прикладных задач, а именно: 

— создание интеллектуальной системы измерения давления для семейства двигате- 
лей АИ-436 самолетов АН-148, Ту-334 [12]; 
— создание ЭСУ температурным режимом двигателя АЙ-25ТЛШ самолета Л-39 [13], 
а также для двигателей семейства ТВЗ-117. 

Перспективы дальнейших исследований заключаются в создании ММ бортово- 

го применения. 
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В.Ф. Миргород, В.М. Грудйиюн 

Нов! комп’ютерно-орентован! форми математичних моделей аващйних газотур@нних двигунв 
Для кола задач оптимального керування, ощнки техн!1чного стану та дагностики автащйних 
газотурбнних двигунв розглянут! питання удосконалення 1снуючих кусково-лйних диференщальних 
та побудови нових 1нтегральних математичних моделей, адекватно в1дображаючих статику та 
динамжу об’ектв керування та д1агностики. Запропонован! нов! математичн! моделг у вигляд! 
1нтегральних равнянь Вольтерра та методи 1х комп’ютерно! реалзаци. 


Г.Е. Мизотов, Г.М. Стиат ип 

'Тве №е Сотрщшег-Онешеа Еогтз$ о? Маетайса! Моде оЁ Амайоп Саз Тигбше Епоше$ 

Еог фе орйта| сопёго| ргоегл$, ап езИтайоп оЁ{есриса| 5е ап4 Ч1а2по$Яс$ оЁ амаНоп газ иифше 
епешез даезНопз оЁ ипргоуетепи оЁ Кпо\уп Ч1Еегепна! ап4 деуеортеп{ ое пе\ п(еогае тафетайса] 
по4е!$ адедиае!у геНесипе айс ап 4упаплс$ оР сопго! ап Фаепто$Нс$ об]ес{ аге ехатшед. Ме\ 
та®фетайса! п1о4е!$ ш бе Ююгт оё пцеотайеа Уо[егга едаайоп$ ап пефо4$ оЁ Феи сотарщег геайтайоп 
аге оНегед. 
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